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Abstract 

The linear stability analysis of the equations governing the evolution of a flat sand 
bed submitted to a turbulent shear flow predicts that the wavelength A at which 
the bed destabilises to form dunes should scale with the drag length Ld ra g = ^ d. 
This scaling law is tested using existing and new measurements performed in water 
(subaqueous ripples), in air (aeolian dunes and fresh snow dunes), in a high pressure 
CO2 wind tunnel reproducing conditions close to the Venus atmosphere and in the 
low pressure CO2 martian atmosphere (martian dunes). A difficulty is to determine 
the diameter of saltating grains on Mars. A first estimate comes from photographs 
of aeolian ripples taken by the rovers Opportunity and Spirit, showing grains whose 
diameters are smaller than on Earth dunes. In addition we calculate the effect of 
cohesion and viscosity on the dynamic and static transport thresholds. It confirms 
that the small grains visualised by the rovers should be grains experiencing saltation. 
Finally, we show that, within error bars, the scaling of A with Ldrag holds over almost 
five decades. We conclude with a discussion on the time scales and velocities at which 
these bed instabilities develop and propagate on Mars. 
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Aeolian sand dunes form appealing and photogenic patterns whose shapes 
have been classified as a function of wind regime and sand supply [1,2]. These 
dunes have martian 'cousins' with very similar features [3]. Grains can also 
be transported by water flows, and bed instabilities are commonly reported in 
flumes, channels and rivers [4,5,6]. The same dynamical mechanisms control 
the formation of aeolian (both on Earth and Mars) dunes and subaqueous 
ripples 1 from a flat sand bed. In a first section we describe the framework 
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in which this instability can be understood. A unique length scale is involved 
in this description, namely the sand flux saturation length L sat , which scales 
on the drag length Ldrag — ^d, where d is the grain diameter and p s /Pf the 
grain to fluid density ratio. It governs the scaling of the wavelength A at which 
dunes or subaqueous ripples nucleate. 

However, attempting to test the scaling laws governing aeolian sand transport 
and dunes characteristic size, one faces a major problem: the grain size does 
not vary much from place to place at the surface of Earth and the grain to air 
density ratio is almost constant. This narrow range of grain size composing 
these dunes has a physical origin: large grains are too heavy to be transported 
by the wind and very small ones - dust - remain suspended in air. As shown by 
Hersen et al. [7,8], small-scale barchan dunes can be produced and controlled 
under water, with a characteristic size divided by 800 - the water to air density 
ratio - with respect to aeolian sand dunes. The martian exploration and in 
particular the discovery of aeolian patterns (ripples, mega-ripples and dunes) 
give the opportunity to get an additional point to check transport scaling laws, 
as the martian atmosphere is typically 60 to 80 times lighter than air. Finally, 
data from high pressure CO2 wind tunnel reproducing conditions close to the 
Venus atmosphere [9] as well as fresh snow barchan pictures [10,11] nicely 
complement and confirm the scaling relation. 

A difficulty already risen in [12] is to determine the diameter of saltating grains 
on Mars. A first estimate is derived from the analysis of photographs of the soil 
taken by the rovers Opportunity and Spirit next to ripple patterns, showing 
grains whose diameters are smaller than on Earth dunes [13]. We derive in 
the appendix the effect of cohesion and viscosity on the dynamic and static 
transport thresholds. It confirms that the small grains visualised by the rovers 
are experiencing saltation. Finally, we show that, within error bars, A scales 
with Ldrag over almost five decades and discuss the corresponding time scales 
and velocities at which these bed instabilities develop and propagate. 



1 Dune instability 

The instability results from the interaction between the sand bed profile, which 
modifies the fluid velocity field, and the flow that modifies in turn the sand 
bed as it transport grains. Let us consider a flat and symmetric bump as 
that shown in figure 1. The fluid is accelerated on the upwind (stoss) side 
and decelerated on the downwind side. This is schematized by the flow lines 
in this figure: they are closer to each other above the bump. This results in 
an increase of the shear stress r applied by the flow on the stoss side of the 

that the water height, whereas that of subaqueous 'dunes' is comparable to depth. 
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Fig. 1. Schematic of the instability mechanism showing the stream lines around a 
bump, the fluid flowing from left to right. A bump grows when the point at which 
the sand flux is maximum is shifted upwind the crest. The shift of the maximum 
shear stress scales on the size of the bump. The spatial lag between the shear and 
flux maxima is the saturation length L sat . 

bump. Conversely, r decreases on the lee side. Assuming that the maximum 
amount of sand that can be transported by a given flow - the saturated sand 
flux - is an increasing function of r, erosion takes place on the stoss slope 
as the flux increases, and sand is deposited on the lee of the bump. If the 
velocity field was symmetric around the bump, the transition between erosion 
and deposition would be exactly at the crest, and this would lead to a pure 
propagation of the bump, without any change in amplitude (we call this the 
'A' effect, see below). In fact, due to the simultaneous effects of inertia and 
dissipation (viscous or turbulent), the velocity field is asymmetric (even on a 
symmetrical bump) and the position of the maximum shear stress is shifted 
upwind the crest of the bump (the '£>' effect). In addition, the sand transport 
reaches its saturated value with a spatial lag, characterized by the saturation 
length L sat . The maximum of the sand flux q is thus shifted downwind the 
point at which r is maximum by a typical distance of the order of L sat . The 
criterion of instability is then geometrically related to the position at which the 
flux is maximum with respect to the top of the bump: an up-shifted position 
leads to a deposition of grains before the crest, so that the bump grows. 

The above qualitative arguments can be translated into a precise mathematical 
framework, see e.g. [14,15]. For a small deformation of the bed profile h(t,x), 
the excess of stress induced by a non-flat profile can be written in Fourier 
space as: 

r = T (A + iB)kh, (1) 
where To is the shear that would apply on a flat bed and k is the wave vector 
associated to the spatial coordinate x. A and B are dimensionless functions 
of all parameters and of k in particular. The A part is in phase with the bed 
profile, whereas the B one is out of phase, so that the modes of h and r of 
wavelength A have a spatial phase difference of the order of XB /(2irA). This 
shift is typically of the order of 10% of the length of the bump as A and B are 
typically of the same order of magnitude. Expressions for A and B have been 
derived by Jackson and Hunt [16] in the case of a turbulent flow. As shown 
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Fig. 2. Sand transport relaxation lengths in units of Ld ra g as a function of the 
rescaled wind velocity u^/u^, as predicted by the theoretical model presented in 
[20]. The inset shows the same curves in linear scale. Starting from a vanishing 
flux, the sand transport first increases exponentially over a distance shown with o 
symbols. The dominant mechanism during this phase is the ejection of new grains 
when saltons collide the sand bed. The spatial corresponding growth rate diverges 
at the threshold and rapidly decreases with larger u*. As the number of grains 
transported increases, the wind is slowed down in the saltation curtain, until satu- 
ration. The distance over which the flux relaxes towards its saturated value is shown 
with • symbols. Except very close to the threshold, the dominant mechanism is the 
negative feedback of transport on the wind. 

by Kroy et al. [17], A and B only weakly (logarithmically) depend on the 
wavelength so that they can be considered as constant for practical purpose. 

If the shear stress is below the dynamical threshold r th , no transport is ob- 
served, hence the sand flux is null. Above this threshold, one observes on a 
flat bed a saturated flux Q which is a function of tq. The fact that a wind 
of a given strength can only transport a finite quantity of sand is due to the 
negative feedback of the grains transported by the fluid on the flow velocity 
profile - the moving grains slow down the fluid. This saturation process of 
the sand flux is still a matter of research [18,19,20]. We refer the interested 
reader to [20,21] for a review discussion. For our present purpose, we need a 
first order but quantitative description of the saturated flux. Following [20], 
Rasmussen et al. wind tunnel data [22] are well described by the relationship: 



g is gravitational acceleration. The prefactor 25 has been adjusted to fit the 
data and is reasonably independent of the grain size d. The equivalent of the 
relation (1) for the saturated flux on a modulated surface g sat can then be 




To > r th and Q — otherwise. 



(2) 
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written as 

q Bat = Q(A + iB)kh, (3) 

where, by the use of relation (2), the values of A and B simply verify A/A = 
B/B = l-r th /r . 

As for any approach to an equilibrium state, there exists a relaxation length 
- or equivalently a relaxation time - scale associated with the sand flux satu- 
ration. This was already mentioned by Bagnold who measured the spatial lag 
needed by the flux to reach its saturated value on a flat sand patch [23]. A 
saturation length L sat in dune models has been first introduced by Sauermann 
et al. [24], where the dependence of L sat on r and in particular its divergence 
as t — > r t h has been put phenomeno logically in the description. In fact, the 
saturation length should a priori depend on the mode of transport - at least 
we are sure that it must exist in all the situations where there is an equilibrium 
transport. As there can be different mechanisms responsible for a lag before 
saturation (delay to accelerate the grain to the fluid velocity, delay due to the 
sedimentation of a transported grain [25], delay due to the progressive ejection 
of grains during collisions, delay for the negative feedback of the transport on 
the wind, delay for electrostatic effects between the transported grains and 
the soil, etc), the dynamics is dominated by the slowest mechanism so that 
L sat is the largest among the different possible relaxation lengths. In the aeo- 
lian turbulent case, there exists a detailed theoretical analysis [20] providing 
the dependence of the saturation length on the shear velocity (• in figure 2). 
The curve roughly presents two zones. Very close to the threshold, the slowest 
process is the ejection of grains during collision. It is thus natural to have a 
divergence of the saturation length at the threshold [24] as the replacement ca- 
pacity crosses 1 (see the calculation of the dynamical threshold in Appendix). 
From just above the threshold - say for u*/u t h > 1.05 - the saturation length 
gently increases (roughly linearly) with u*. However, in the field, the mean 
wind strength varies from day to day, as seasons goes by. For practical pur- 
poses - e.g. photograph analysis - it is thus of fundamental importance to 
define an effective saturation length, independent of the wind velocity. Fortu- 
nately, the velocity rarely exceeds ~ 3 u t h, so that, in the range of interest, the 
velocity dependence on L sat is subdominant. Our first important conclusion is 
thus that even though it should be remembered that L sat slightly depends on 
u^/u t h in lab experiments, the dominant parameters are the grain size d and 
the sand to fluid density ratio. The theoretical analysis [20] shows that L sat 
scales on the drag length Ldrag = ^ d, which is the length needed for a grain 
in saltation to be accelerated to the fluid velocity. It is worth noting that this 
inertial effect is however not the mechanism limiting the saturation process 
(see above). Using the results of a field experiment performed in the Atlantic 
Sahara of Morocco [15], the prefactor between L sat and L drag can be computed 
and gives: 

L sat ~4.4^d. (4) 
Pf 



5 



W 1D ' 



io 2 S (unn) 




Fig. 3. Measurement of the typical grain size from the auto-correlation function 
C(S) of a photograph of the granular bed. a Sample of aeolian sand from the 
mega-barchan of Sidi-Aghfinir (Atlantic Sahara, Morocco), of typical diameter 
d — 165 /im. Top inset: Probability distribution function of the grain diameter 
d, weighted in mass. The narrow distribution around the maximum of probability is 
characteristic of the aeolian sieving process. Bottom inset: 200 x 200 pix 2 zoom on 
a typical photograph on which the computation of the autocorrelation function has 
been performed. The image resolution is here 5 pix/d. b Sample of sand from Mars 
(at rover Spirit's landing site: 16.6° S, 184.5° W). Top inset: photograph showing 
the presence of milimetric grains as well as much smaller ones which are expected to 
be the saltons. Bottom inset: same as a, but the typical image resolution is 3 pix/d. 



This result is also consistent with Bagnold's data [23]. Note that the saturation 
length has never been measured directly in other situations (neither under 
water nor in high/low pressure wind tunnels). 

The linear stability analysis of the coupled differential equations of this frame- 
work has been performed in [14] . In particular, the wavenumber corresponding 
to the maximum growth rate is given by 



1rt X 1 / 3 3B 

Aw^sat = X 1/3 — with X = - ir 

3 A 



3 + J3(l + (A/By) 



(5) 



For typical values of A and B determined in the barchan dune context [15], 
we have A max = 27r/k max ~ 12L sat . Note that we have A/B = A/B, so that 
/c max L sat is independent of r and r th . Using measurements in water (subaque- 
ous ripples), in air (aeolian dunes) and in the CO2 atmosphere of Mars and of 
the Venus wind tunnel, we now investigate the scaling relation between A max 
and L drag . To do so, we need to first solve the controversy concerning the size 
of saltating grains on Mars. 
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2 Size of saltating grains on Mars 

The determination of the typical size of the grains participating to saltation on 
martian dunes is a challenging issue. First, the dunes seem to evolve very slowly 
or may even have become completely static. Second, no sample of the matter 
composing the bulk of the dunes is available. Third, we do know from the 
observation of dunes on Earth that they can be covered by larger grains that do 
not participate to the transport in saltation. With the two rovers Opportunity 
and Spirit, we now have direct visualisations of the soil [13], and in particular 
of clear aeolian structures like ripples 2 or nabkhas 3 . Unfortunately, these 
structures did not lie on the surface of a dune. We thus analyze the available 
photographs, assuming that, like on Earth, the size of the grains participating 
in saltation does not vary much from place to place. 



2. 1 Direct measurement of grain sizes 

The photographs freely accessible online are not of sufficiently good resolu- 
tion to determine the shape and the size of each grain composing the surface. 
Besides, one has to be careful when analyzing such pictures as part of what is 
visible at the surface corresponds to grains just below it and partly hidden by 
their neighbours. We have thus specifically developed a method to determine 
the average grain size in zones where the grains are reasonably monodispersed, 
when the resolution is typically larger than 3 pixels per grain diameter. This 
measure can be deduced from the computation of the auto- correlation func- 
tion C{$) of the picture, which decreases typically over one grain size. More 
precisely, we proceed with the following procedure: 

• The zones covered by anything but sand (e.g. gravels or isolated larger 
grains) are localized and excluded from the analysis. 

• Because in natural conditions the light is generally inhomogeneous, we per- 
form a local smoothing of the picture with a gaussian kernel of radius ~ 10 d. 
The resulting picture is subtracted from the initial one. 

• We compute the local standard deviation of this image difference with the 
same gaussian kernel and produce, after normalisation, a third picture X of 
null local average and of standard deviation unity. 

• The auto-correlation function C(S) is computed on this resulting picture, 



2 Contrarily to dunes, aeolian sand ripples form by a screening instability related 
to geometrical effects [23,26]. 

3 As there is a reduction of pressure in the lee of any obstacle, sand accumulates in 
the form of streaks aligned in the direction opposite to the wind (shadow dunes). 
These structures are called nabkhas [27]. 
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averaging over all directions: 



J2k,l/k 2 +l 2 =5 2 1 



Figure 3a shows the curve C{5) obtained from a series of photographs of aeo- 
lian sand sampled on a terrestrial dune. For all the resolutions used (between 
1 and 5 pix/d), the data collapse on a single curve, which is thus characteristic 
of the sand sample. The top inset shows the distribution of size, weighted in 
mass, in log-log representation. It presents a narrow peak around the value. 
The autocorrelation function decreases from 1 to over a size of the order of 
d 50 (Lorenzian fit). This is basically due to the fact that the color or the gray 
level at two points are only correlated if they are inside the same grain. It is 
thus reasonable to assume that the autocorrelation curve is a function of 5/d 5Q 
only. As a matter of fact, photographs of two samples of comparable polydis- 
persity are similar once rescaled by the mean grain diameter. We will thus use 
the curve obtained with aeolian sand sampled on Earth as a reference to de- 
termine the size of Martian grains. Figure 3b shows the autocorrelation curve 
obtained from the colour image taken by the rover Spirit at its landing site. 
C(S) decreases faster than the reference curve but by tuning the value of the 
martian d 50 , one can superimpose the Mars data with the solid line fairly well. 
From this picture as well as a series of gray level photographs taken by the 
rover Opportunity, we have estimated the diameter of the grains composing 
these aeolian formations to 87 ± 25 jum. 

The first measurement of martian grain sizes dates back to the Viking mission 
in the 70s. On the basis of thermal diffusion coefficient measurements, Edgett 
and Christensen have estimated the grain diameter to be around 500 /mi and 
at least larger than those composing dunes on Earth [12]. This size corresponds 
to larger grains than saltons such as those shown in the top inset of figure 3b. In 
agreement with our findings, wind tunnel experiments performed in 'martian 
conditions' [28] have shown that grains around 100 fim are the easiest to 
dislodge from the bed. We shall come back to this point later on. 

Several theoretical investigations of the size of aeolian martian grains have 
been conducted, starting with the work of Sagan and Bagnold [29]. In that 
paper the authors argue that, as Mars is a very arid planet, the cohesion forces 
due to humidity can be neglected and proposed a cohesion-free computation 
which predicts that very small grains (typically of one micron) may be put into 
saltation. Miller and Komar [30] also followed this cohesion-free approach and 
proposed static threshold curves with no turnup on small particle size. It was 
soon realized that cohesion forces can occur for reasons other than humidity - 
namely van der Waals forces - and several authors proposed (static) transport 
threshold curves with a peaked minimum around 100 fim [31,32,33,34,35]. 
However in these papers, cohesion is treated in an empirical fashion, with the 
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Fig. 4. Diagram showing the mode of transport on Mars as a function of the grain 
diameter d and of the turbulent shear velocity « t h (left) or of the wind speed Uth at 
2 m above the soil (right). Below the dynamical threshold (dashed line), no grain 
motion is observed. A grain at rest on the surface of the bed starts moving, dragged 
by the wind, when the velocity is above the static threshold (solid line). Between the 
dynamical and static thresholds, there is a zone of hysteresis where transport can 
sustain due to collision induced ejections. The progressive transition from saltation 
to suspension as wind fluctuations become more and more important is indicated 
by the gradient from white to gray color. See appendix for more details on the 
derivation of this graph. 

assumption that van der Waals forces lead to an attractive force proportional 
to the grain diameter with a prefactor independent of d. 

We have recomputed the Martian transport diagram (figure 4) using a new 
derivation of the transport thresholds. It takes into account the hysteresis 
between the static and dynamic thresholds, the effect of viscosity and - in 
a more rigorous way - the effect of cohesion. As this derivation, although 
consequent and original, is not the central purpose of the present paper which 
is devoted to the test of the dune scaling law, we have developed and discussed 
it in appendix. 

We wish to solely discuss here the figure 4, which is useful to prove that the 
87 fim sized grains can be transported in saltation. The first striking feature 
of the Martian transport diagram is the huge hysteresis between the dynamic 
and static thresholds. Compared to aeolian transport on Earth (figure A. 2 
in the Appendix) for which the static threshold is typically 50% above the 
dynamic threshold, they are separated on Mars by a factor larger than 3. 
Quantitatively, if the static threshold is very high compared to the typical 
wind velocities on Mars (~ 150 km/h at 2 m above the soil), the dynamical 
threshold is only of the order of ~ 45 km/h at 2 m. From images by the 
Mars Orbiter Camera (MOC) of reversing dust streaks, Jerolmack et al. [36] 
have estimated that modern surface winds can reach velocities as large as 
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Table 1 

Comparison of different quantities (gravity g, initial wavelength of bed instability 
A, diameter of saltons d as well as fluid and sediment densities pf and p s ) in the 
air (Earth), in the martian and Venus wind tunnel CO2 atmospheres and in water. 
As the temperature at the surface of Mars can vary by an amplitude of typically 
100 K between warm days and cold nights, the density of the atmosphere displays 
some variation range. 

150 km/h. Looking at figure 4, it can be seen that at such velocities, most of the 
grain below 100 /zm can be suspended and that even millimetric grains can be 
entrained into saltation. Even in less stormy conditions, sand transport should 
not be as unfrequent as one could expect, even though the large amplitude of 
the hysteresis implies an intermittency of sand transport. This suggests that 
Martian dunes are definitively active. 



3 A dune wavelength scaling law 

Keeping in mind that the wavelength A that spontaneously appears when a 
flat sand bed is destabilized by a turbulent flow scales on the flux saturation 
length, the aim of the paper is to plot A against Ld rag in the different situations 
mentioned in the introduction: aqueous ripples, waves on aeolian dunes, both 
on Earth and Mars, fresh snow dunes and microdunes in the Venus wind 
tunnel - numerical values of the parameters corresponding to these different 
situations are summarized in table 1. 

Our reference data point in figure 5 (aeolian dunes on Earth) has been ob- 
tained after an extensive work on barchan dunes in the Atlantic Sahara of 
Morocco. In a recent paper [15], we have shown that perturbations such as 
wind changes generate waves at the surface of dunes by the linear instabil- 
ity described above. Direct measurements of the wavelength are reported in 
Figure 6a in an histogram. They give A ~ 20 m for waves on the flanks of 
medium sized dunes (solid line) and ~ 28 m on the windward side of a mega- 
barchan (dashed line), where some pattern coarsening occurs. L sat ~ 1.7 m 
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Fig. 5. Average wavelength as a function of the drag length. The two black squares 
(resp. diamonds) labeled 'Earth' (6) (resp. 'Mars' (o*)) come from the two different 
histograms of figure 6a (resp. 6b). The three white circles are the under water 
Coleman and Melville's data [5,6], whereas the black circle is that from Hersen et 
a/.'s experiments [7,8]. The white triangle has been computed from the Venus (?) 
wind tunnel study [9]. Snow (-X-) dune photos (see figure 8) have been calibrated 
and complemented with data from [10,11,37], and give the white squares. 
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Fig. 6. Histograms of the wavelength measured in the Atlantic Sahara (Morocco) 
and in Mars southern hemisphere (mainly, but not only, in the region 320-350° W, 
45-55° S). (a) Earth (s). Solid line: wavelengths systematically measured on barchan 
dunes located in a 20 km x 8 km zone; Dash line: wavelengths measured on the 
windward side of a mega-barchan whose surface is permanently corrugated and 
where some coarsening occurs, (b) Mars ((f). Solid line: wavelengths measured on 
dunes in several craters (e.g. Rabe, Russell, Kaiser, Proctor, Hellespontus); Dash 
line: histogram restricted to Kaiser crater (341°W, 47°S). Averaged values of A are 
20 m (solid line on panel a), 28 m (dash line on panel a), 510 m (solid line on panel 
b) and 606 m (dash line on panel b). 
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Fig. 7. Comparison of dune morphologies on Earth and on Mars, a Mega-barchan 
in Atlantic Sahara, Morroco. b 'Kaiser' crater on Mars (341°W,47°S). c Closer view 
of the Kaiser crater dunes. As on Earth, small barchans are visible on the side of 
the field. 




Fig. 8. Fresh snow aeolian dunes on ice. (a) Transverse dunes on the iced baltic 
sea (credits Berth Hakansson, Swedish Meteorological and Hydrographical Insti- 
tute/Baltic Air-Sea-Ice Study). The wavelength is around 15 m [37]. (b) Snow 
barchan field in Antartica (credits Stephen Price, Department of Earth and Space 
Sciences, University of Washington) . The shadow of the Twin Otters, which has has 
a wing span of 19.8 m and a length of 15.7 m, gives the scale. In both pictures, the 
perspective has been corrected to produce a pseudo-aerial view. 
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was measured independently in the same dune field, and these aeolian grains 
of 180 (Ltm lead to a drag length of 40 cm. 

The study of small scale barchans under water [7,8] provides the second data 
point: measurements give A ~ 2 cm with glass beads of size d = 150 /zm, a 
size which leads to a drag length of 400 fim. This point is the only one in the 
literature for which it is specified that: (i) the wavelength has been measured 
during the linear stage of the instability; (ii) the sand bed destabilizes in 
a homogeneous way and not starting at the entrance of the set-up, due to 
strong disturbances there; (iii) the height of the tank is much larger than the 
wavelength. We have added three other underwater points in figure 5, from 
Coleman and Melville's data [5, 6]. The slight discrepancy with the straight 
line (these points are all above the line) can be explained by the three points 
risen above: it may be induced from [5] that (i) the initial stage is not clearly 
resolved; (ii) the waves seem to nucleate on defects as if the transition was 
sub-critical; (iii) the flume is only 0.28 m deep, a height which is comparable 
to the observed wavelengths. 

Recent photos of martian dunes [3] such as that in figure 7 lead to an estimate 
of the value of A ~ 600 m on Mars. We focused on the dunes found in craters 
of the southern hemisphere. For comparison, we have measured wavelengths 
on dunes in several craters (e.g. Rabe, Russell, Kaiser, Proctor, Hellespontus) 
and also produced a histogram restricted to Kaiser crater, see figure 6b. As 
for the drag length, we used the value for the grain diameter estimated in the 
previous section. The density of martian grains is similar to or slightly higher 
[36] than that of terrestrial grains. The value of the density of the martian 
atmosphere however varies within some range as the temperature on Mars 
surface can change by an amplitude of typically 100 K between warm days 
and cold nights. In the end, it gives a value of L drag between 13 and 17 m. 

The fourth data point is obtained from Greeeley et al.'s experiment in a high 
pressure C02 wind tunnel [9], which gives for Ldrag a value we expect on 
Venus. The observed wavelength first decreases from 18 cm to 8 cm with the 
wind speed from u* = u t h to u* = 1.6 -u t h and then increases up to 27 cm at 
u* — 2.1 u th . Above this value, the flat sand bed was again found to be stable. 
Although some of these features are reminiscent of the L sat curve discussed 
above (figure 2), this series of experiments is also questionable: nothing is 
reported about the nature of the destabilization (homogeneous appearance of 
the pattern or not) nor on the maturity of the pattern when the wavelength 
was measured (coarsening?). Remarkably, the crude averaged of the provided 
data lies on the master curve, even though the smallest wavelength is less than 
half the predicted value. 

As for the Venus wind tunnel, the fifth data points are not very precise and 
correspond to fresh snow dunes formed on Antartic sea ice and on Baltic sea 
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ice. The typical wavelength A of transverse dunes (for instance in figure 8(a) 
ranges from 15 m to 25 m [37]. Whenever barchan dunes form, they are typi- 
cally 5 m to 10 m long and are separated by ~ 20 m (see figure 8(b). This is 
very similar with aeolian ones, although snow barchan dunes look more elon- 
gated. The determination of Ldrag is more problematic as the snow density 
rapidly increases once fallen. Snow dunes are rather rare and probably form 
only at the surface of ice by strong wind, when the snow can remain fresh and 
not very cohesive. To get the numbers, we have used the measurements per- 
formed close to the Antartic dunes by Massom et al. [10,11]. Although one may 
find somehow disappointing that on the graph of figure 5 these data points are 
located very close to the aeolian sand dunes, they are particularly interesting 
as they show that one can keep the same A by changing simultaneously p s and 
d in a compensating way. 



Globally, we obtain a consistent scaling law A ~ 53Ld rag that covers almost 
five decades (figure 5). We wish to emphasize again that we do not claim 
to capture all the dependences on this plot, but that Ldrag is the dominant 
scaling factor. As shown in figure 2, we expect subdominant dependences on 
the wind speed, on finite size effects, etc... Note that this analysis explains for 
example why martian dunes whose sizes are of the order of a kilometer are 
not 'complex' or 'compound' as their equivalents on Earth [38]. Following our 
scaling relation, they correspond to the small terrestrial dunes and the whole 
dune field on the floor of a crater should rather be considered as a complex 
martian dune. A 'similarity law' for the size of (developed) dunes with L sat 
was schematically drawn by Kroy et al. [39], supported by the existence of 
centimetric barchans under water [7] . However, this similarity was announced 
to fail with Martian dunes. In fact, the scaling of large dunes (i.e. whose sizes 
are much larger than the elementary length A) with L^ Tag (or L sat ) is far from 
obvious as the size selection of barchans involves secondary instabilities related 
to collisions and fluctuations of wind direction [15,40,41]. 



Finally, we would like to end this section with the prediction of the wavelength 
at which a flat sand bed should destabilize on Titan where (longitudinal or 
seif) dunes have been recently discovered [42]. The atmosphere there is ap- 
proximately four times denser than on Earth, and the grains are believed to 
be made of water ice [43]. Computing the threshold curves (not shown), one 
observes that the dynamical and static thresholds are almost identical and 
that the minimum of threshold shear stress is reached for 160 /xm. Using this 
size for the grain diameter a density ratio of the order of 200, we find a cen- 
timetric drag length. Following the scaling law of figure 5, this would lead to 
a dune wavelength between 1 and 2 meters. Note that, unfortunately, this is 
well below the resolution of Cassini radar. 
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4 Discussion: time scales and velocity for the martian dunes 



Now that we have this unique length scale at hand for dunes, it is interesting to 
address the question of the corresponding growth time scales and propagation 
velocities. The linear stability analysis [14,15] tells that the growth rate a 
and the propagation velocity c of bedforms are related to the wavenumber 
k = 2n/\ as 



2 B - AkL sat _ A + BkL sat 

' i + (kL sat y and c{k) - Q i + (kL sat r 



where Q is the saturated sand flux - assumed constant - over a flat bed. 
In order to make these relations quantitative for the bedforms on Earth and 
Mars, we need an effective time averaged flux Q. Recall from equation (2) 
that Q can be related to the shear velocity u*. For simplicity, we suppose that 
there is sand transport (w* > w t h) a fraction r\ of the time and that the shear 
stress is then constant and equal to (1 + a)pfU^ h . The time averaged sand flux 
can thus be effectively expressed as 

g = 25a^/?^<. (8) 
V 9 Ps 

The values of the coefficients A and B which come into equation (7) depend 
on the excess of shear above the threshold as A/ A = Bj B = (1 + a) /a. The 
timescale over which an instability develops is that of the most unstable mode. 
This means that a and c should be evaluated at k — k max (see equation (5)) 
and thus scale as: 



c ^ (i+Q Mtf^ (io> 

Using meteorological data and measurements of the average dune velocities, 
we estimated that Q& is between 60 to 90 m 2 /year [15] and 775 between 65% and 
85%. This gives an effective value aj between 1.5 and 2. Calculating explicitly 
the prefactors, we get a growth time a^ 1 ~ 2 weeks and cj ~ 200 m/year, 
values which are consistent with direct observation. Note that time scales for 
individual fully developed dunes depend on their size, but are also proportional 
to 1/Q [40]. 

In order to compare these terrestrial values to those on Mars, we need to 
estimate the different factors which come into the expressions (9) and (10). 
Recall that, on Earth, for grains of 165 jum, we have U(, th ~ 0.2 m/s. The 
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corresponding value for the 87 /zm sized Martian grains can be found on fig- 
ure 4 and reads Worth ~ 0-64 m/s. The Mars to Earth ratios for the fluid and 
grain densities, as well as the gravity acceleration are known. As discussed in 
the calculation of the transport thresholds (see Appendix) one expects that 
«cf is of order unity and for simplicity we simply take = a&. Finally, the 
most speculative part naturally concerns the value of i]^. Assuming that the 
winds on Mars are similar to those on Earth, we computed on our wind data 
from Atlantic Sahara the fraction of time during which u* > 0.64 m/s and 
got ~ 3%. This value corresponds to few days per year and is probably re- 
alistic as, besides, the soil of Mars is frozen during the winter season. With 
these numerical values, we find that a& is smaller than as by more than five 
decades. In other words, the typical time over which we could see a significant 
evolution of the bedforms on martian dunes is of the order of tens to hundreds 
centuries. Similarly, the ratio Co»/q, is of the order of 10~ 4 . Note that these 
Mars to Earth time and velocity ratios are proportional to i]^, so that larger or 
smaller values for this speculative parameter do not change dramatically the 
conclusions, the main contribution being the density ratio to the power 2 or 
3 (see equations (9-10)). Therefore, as some satellite high resolution pictures 
definitively show some evidence of aeolian activity - e.g. avalanche outlines - 
it may well be that the martian dunes are fully active but not significantly at 
the the human scale. 



5 Conclusion 



As for the last section, we would like to conclude the paper with a summary of 
the status of the numerous hypothesis and facts we have mixed and discussed. 
Although coming from recent theoretical works, it is now well accepted that 
the wavelength A at which dunes form from a flat sand bed is governed by 
the so-called saturation length L sat . However, the dependences of L sat with 
the numerous control parameters (Shields number, grain and flow Reynolds 
numbers, Galileo number, grain to fluid density ratio, finite size effects, etc) 
is still a matter of debate. In this paper, we have collected measurements of 
A in various situations that where previously thought of as disconnected: sub- 
aqueous ripples, micro-dunes in Venus wind tunnel, fresh snow dunes, aeolian 
dunes on Earth and dunes on Mars. We show that the averaged wavelength 
(and thus L sat ) is proportional to the grain size times the grain to fluid den- 
sity ratio. This does not preclude sub-dominant dependencies of L sat with the 
other dimensionless parameters. In each situation listed above, we have faced 
specific difficulties. 

Subaqueous ripples — The transport mechanism under water, namely the di- 
rect entrainment by the fluid, is different from the four other situations. The 
saturation process could thus be very different in this case. The longest re- 
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laxation length could for instance be controlled by the grain sedimentation, 
as recently suggested by Charru [25]. As there is also negative feedback of 
transport on the flow, the destabilization wavelength would be larger than 
our prediction. The formation of subaqueous ripples has remained controver- 
sial mostly because the experimental data are very dispersed. Most of them 
suffer from finite size effects (water depth or friction on lateral walls), from un- 
controlled entrance conditions leading to an inhomogeneous destabilization, or 
from a late determination of the wavelength after a pattern coarsening period. 
For the furthest left data point only, we are certain that all these problems 
were avoided. 

'Venus' micro-dunes — This beautiful experiment gives the most dispersed 
data of the plot. Part of it may be due to the defects listed just above. The com- 
plex variation of the wavelength with the wind speed observed experimentally 
could be interpreted as a sub-dominant dependance of the saturation length 
on the Shields number. However, not only the size of the dunes changes but 
also their shapes. This is definitely the signature of a second length scale at 
work. 

Fresh snow dunes — On the basis of existing photographs, we have clearly 
identified snow dunes pattern resembling sand aeolian destabilization ones. 
They form under strong wind, on icy substrate. The obvious difficulty is that 
the precise state of the snow flakes during the dune formation involves com- 
plex thermodynamical processes that do not exist in the other cases. 
Aeolian dunes on Earth — Although the instability takes place in the field 
under varying wind conditions, the resulting wavelength is very robust from 
place to place and time to time. We consider this point as the reference one 
in the plot. 

Martian dunes — The specific difficulty of the Martian case does not come 
from the wavelength measurements, as the available photographs are well re- 
solved in space, but rather from the prior determination of the diameter d of 
the grains involved. The controversy comes from early determinations based 
on thermal diffusivity measurements, concluding that dunes are covered by 
large (500 /im) grains. Such large grains would lead to a significant shift of 
the data point towards the right of figure 5. A large part of this paper is 
consequently devoted to an independent determination of d, based on the 
analysis of the Martian rovers photographs. Our determination is very differ- 
ent: 87 ± 25 /im. We have computed the sand transport phase diagram with 
more subtleties than in the available literature (including in particular hys- 
teresis and cohesion) and shown that this size corresponds to saltating grains. 
A directly related controversial point was the state of the Martian dunes (still 
active or fossile). We have shown that moderately large winds can transport 
the grains (the dynamical threshold is much lower than the static one) but that 
the characteristic time scale over which dunes form is five orders of magnitude 
larger than on Earth. 



17 



The Moroccan wavelength histogram has been obtained in collaboration with 
Hicham Elbelrhiti. We thank Eric Clement and Evelyne Kolb for useful ad- 
vices on the way cohesion forces can be estimated. We thank Frangois Charru 
and Douglas Jerolmack for discussions. We thank Brad Murray for a careful 
reading of the manuscript. Part of this work is based on lectures given during 
the granular session of the Institut Henri Poincare in 2005. This study was 
financially supported by an 'ACI Jeunes Chercheurs' of the french ministry of 
research. 



18 



A A model for transport thresholds, including cohesion 

We have directly measured the grain diameter (87 ±25 /im) in Martian ripples 
photographed by the rovers. We then made the double assumption that (i) 
the grains composing the martian dunes are of the same size and that (ii) 
these grains participate to saltation transport whenever the wind is sufficiently 
strong. In order to support hypothesis (ii), we have computed the transport 
thresholds in the Martian conditions. We find that the grains which move first 
for an increasing wind are around 65 /im in diameter. The full discussion of the 
transport thresholds is however too heavy to be incorporated into the body 
of the article, mainly devoted to the dune wavelength scaling law. We give 
below a short but self-sufficient derivation of the scaling laws for the static 
and dynamic transport thresholds. 

A.l Definition of transport thresholds, hysteresis 

We consider the generic case of a fluid boundary layer over a flat bed composed 
of identical sand grains. For given grains and surrounding fluid, the shear stress 
r controls the sand transport. The dominant mechanism for grain erosion de- 
pends on the sand to fluid density ratio. In dense fluids grains are directly 
entrained by the flow whereas in low density fluids grains are mostly splashed 
up by other grains impacting the sand bed. Two thresholds are associated to 
these two mechanisms: starting from a purely static sand bed, the first grain 
is dragged from the bed and brought into motion at the static threshold r sta . 
Once sand transport is established, it can sustain by the collision/ejection 
processes down to a second threshold Td yn . The sand transport thus presents 
an hysteresis - responsible for instance for the formation of streamers. Fi- 
nally, when the fluid velocity becomes sufficiently high, more and more grains 
remain in suspension in the flow, trapped in large velocity fluctuations. Al- 
though there is no precise threshold associated to suspension, one expect that 
these fluctuations becomes dominant for the grain trajectories when the shear 
velocity u* = \Jr/p~f is much larger than the grain sedimentation velocity u^n. 
The behaviour of these thresholds in the Martian atmosphere, as well as in 
water and air, with respect to the grain diameter is summarized in figures 4 
and A. 2, and we shall now discuss how to compute these curves. 

A. 2 An analytic expression for the turbulent boundary layer velocity profile 

In the dune context, the flow is generically turbulent far from the sand bed 
and the velocity profile is known to be logarithmic. However, our problem is 
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more complicated as we have to relate the velocity u of the flow around the 
sand grains (i.e. the velocity at the altitude, say, z = d/2 that we call v in the 
following) to the shear velocity u*. Depending on the grain based Reynolds 
number, there either exists a viscous sub-layer between the soil and the fully 
turbulent zone, or the momentum transfer is directly due to the fluctuations 
induced by the soil roughness. The first step is thus to derive an expression 
for the turbulent boundary layer wind profile valid in the two regimes, far or 
close to the bed. To do so, we express the shear stress as the sum of a visous 
and a turbulent part as 



r = pfi>d z u + k 2 (z + rd) 2 pf\d z u\d z u, 



(A.I) 



where v is the kinematic viscosity of the fluid, k the von Karman constant 
and r the aerodynamic roughness rescaled by the grain diameter d. The shear 
stress is constant all through the turbulent boundary layer and equal to PfU 2 . 
Let us define the grain Reynolds number as: 



Ren = 



2Kii*rd 



v 



(A.2) 



and a non-dimensional distance to the bed: 



rd 



(A.3) 



Note that Z is tends to 1 on the sand bed and that it is equal to 1 + l/2r 
at the center of the grain. With this notations, the differential equation (A.l) 
can be rewritten under the form: 



d(Ku/u*) 
dZ 



+ 



Ren 



d(nu/u*) 
dZ 



1 = 0, 



(A.4) 



which is easily integrated into: 



u 



K 



sink' 1 (Re Z) + 



1-Jl + RelZ* 



RenZ 



z=z 



Z=l 



(A.5) 



Performing expansions in the purely viscous and turbulent regimes, one can 
show that a good approximation of the relation between the shear velocity u* 
and typical velocity around the grain v is given by 



2u 



K 



U * d V+ ln 2 (l + l/2r) 



v 2 . 



(A.6) 
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Fig. A.l. a Grain packing geometry considered for the computation of the static 
transport thresholds, b Schematic drawing showing the contact between grains at 
the micron scale. 

A. 3 Static threshold: influence of Reynolds number 



The bed load (tractation) threshold is directly related to the fact that the 
grains are trapped in the potential holes created by the grains at the sand bed 
surface. To get the scaling laws, the simplest geometry to consider is a single 
spherical grain jammed between the two neighbouring (fixed) grains below it, 
see figure A. la. Let us first discuss the situation in which the cohesive forces 
between the grains are negligible and the friction at the contacts is sufficient to 
prevent sliding. It can be inferred from figure A. la that the loss of equilibrium 
occurs for a value of the driving force F proportional to the submerged weight 
of the grain: F oc (p s — pf)gd 3 , where p s is the mass density of the sand grain, 
and pf is that of the fluid - which is negligible with respect to p s in the case 
of aeolian transport. As F is proportional to the shear force rd 2 exerted by 
the fluid, the non-dimensional parameter controlling the onset of motion is 
the Shields number, which characterises the ratio of the driving shear stress 
to the normal stress: 

(Ps-Pf)gd' A ' 7 

The threshold value can be estimated from the geometry of the piling, and 
depends on whether rolling or lifting is the mechanism which makes the grain 
move. Finally, the local slope of the bed modifies the threshold value as traps 
between the grains are less deep when the bed is inclined. In particular, its 
value must vanish as the bed slope approaches the (tangent of the) avalanche 
angle. We here ignore these refinements which can be incorporated into the 
values of A and B (see section 1). 

At the threshold, the horizontal force balance on a grain of the bed reads 

\p{p s ~ Pf)gd 3 = \c dPf v 2 th d\ (A.8) 

where p, is a friction coefficient, and Cd the drag coefficient which is a function 
of the grain Reynolds number. With a good accuracy, the drag law for natural 
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grains can be put under the form: 
with Coo — 1 and s ~ 5 for natural sand grains [44] . 



(A.9) 



At this stage, we introduce the viscous size d u , defined as the diameter of 
grains whose free fall Reynolds number u^wd/u is unity: 



du = (Ps/Pf - ir 1/3 ^ g- 1/3 



(A.10) 



It corresponds to a grain size at which viscous and gravity effects are of the 
same order of magnitude. We define v as the fluid velocity at the scale of the 
grain normalised by (p s /pf — l) 1 / 2 ^) 1 / 2 . From the three previous relations, 
we get the equation for v t h' 



1/2 



o, 



(A.ll) 



which solves into 



Vth 



4C. 



3/2 



d 



+ 



fxC 



1/2 X 



1/2 



du_ 
d 



3/4 



(A.12) 



The expression of the static threshold Shields number is finally: 



3/2 



5th + 



ln 2 (l + l/2r) 



v 2 



(A.13) 



A. 4 Static threshold: influence of cohesion 



For small grains, the cohesion of the material strongly increases the static 
threshold shear stress. Evaluating the adhesion force between grains is a dif- 
ficult problem in itself. We consider here two grains at the limit of separation 
and we assume that the multi-contact surface between the grains has been 
created with a maximum normal load N max , see figure A. lb. The adhesion 
force A^h can be expressed as an effective surface tension 7 times the radius 
of curvature of the contact, which is assumed to scale on the grain diameter: 

iVadhOCTU (A. 14) 
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This effective surface tension is much smaller than the actual one, 7, as the 
real area of contact A real is much smaller than the apparent one An eTtz .: 



7 = 7^- (A.15) 

^iHcrtz 

The apparent area of contact can be computed following Hertz law for two 
spheres in contact under a load iV max : 

A Hertz oc (^^.y 3 , (A. 16) 

where E is the Young modulus of the grain [45]. To express the real area of 
contact, we need to know whether the micro- contacts are in an elastic or a 
plastic regime. Within a good approximation, A Tea _\ can expressed in both cases 
[46] as: 

^eal « (A.17) 

where M is the Young modulus E (elastic regime) or the hardness H (plastic 
regime) of the material. Altogether, we then get: 

£2/3 

iVadhOC 7 ^-(iV m axrf) 1/3 . (A.18) 



In order to bring into motion such grains, the shear must be large enough to 
overcome both weight and adhesion, so that, for iV max ~ p s gd 3 (i.e. the weight 
of one grain) the critical Shields number is the sum of two terms and takes 
the form of 



6 



with 



tii 




oc 



7 \ 3 / 5 / Hi \ 



E 



\ 2 /5 



(A.19) 



(A.20) 



M) \p s g) 

Note that, by contrast to the references [31,32,33,34,35], we find that the 
adhesive force finally scales, for d — > with d 4 / 3 and the critical Shields 
number with d~ 5 ^ 3 (instead of exponents 1 and —2 respectively). 



A. 5 Dynamical threshold 



The dynamical threshold can be defined as the value of the control parameter 
for which the saturated flux vanishes. Fitting the flux vs shear velocity relation, 
one can measure the dynamical threshold in a very precise way (see for instance 
the analysis in [20] of the data obtained by Rasmussen et al. [22]). 
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Fig. A. 2. Diagram showing the mode of transport in the aeolian (a) and underwater 
(b) cases, as a function of the grain diameter d and of the turbulent shear velocity 
u t h (left) or of the wind speed U t h at 2 m above the soil (right). The dynamical 
threshold (dashed line) is below the static threshold (solid line) in the aeolian case 
but much above it underwater. The dark gray is the zone where no transport is 
possible. Above, the background color codes for the ratio u*/uf a \i: white corresponds 
to negligible fluctuations and gray to suspension. The experimental points are taken 
from (O) Chepil [48] and (□) Rasmussen [22,20] in the aeolian case and from (O) 
Yalin and Karahan [49] in the underwater case. 

As shown in [20] , the wind velocity profile is almost undisturbed at the dynam- 
ical threshold. The shear stress threshold is thus achieved when the velocity 
that a grain acquires after a single jump is just sufficient to eject on average 
one single grain out of the bed after a collision (unit replacement capacity cri- 
terion). The impact velocity at this dynamical threshold is thus proportional 
to the trapping velocity i.e. the velocity needed for a grain to escape from its 
potential trapping (see Quartier et al. [47]): 



An analytical expression of the dynamical threshold has been derived in [20], 
but only in the limit of large Reynolds numbers. To generate the plots dis- 
played in figures 4 and A. 2, we use here a numerical integration of the equa- 
tions of motion of a grain - with equations (A. 5) and (A. 9) for the wind profile 
and drag coefficient. The trajectories are computed for an ejection velocity v-\ 
equal to evi and a typical ejection angle 7r/4 (see [20]) . 

A. 6 From Earth to Mars 

In order to get a reasonable transport diagram on Mars, we have first tuned the 
different parameters controlling the dynamic and static thresholds to repro- 




(A.21) 
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duce within the same model the underwater and aeolian data (figure A. 2). For 
the friction coefficient fi, we have simply taken the avalanche slope for typical 
aeolian grains tan(32°). As already found in [20], the rescaled soil roughness r 
is higher than the value found by Bagnold. It is here adjusted to r = 1/10. The 
value of a, the impact velocity needed to eject statistically one grain rescaled 
by the trapping velocity, is adjusted to 15. The restitution coefficient was ad- 
justed to e = 2/3. For obvious reasons it is possible to obtain almost the same 
curves with different pairs (a, e). Finally the cohesion length d m was tuned to 
25 /an, which is consistent with the above calculation. 

We see in figure A. 2 that the agreement is excellent both in the aeolian and 
underwater cases. Due to the low density ratio in water, the collision process 
is completely inefficient. As a consequence, the dynamical threshold is well 
above the static one: the only erosion mechanism is the direct entrainment 
of grains by the fluid (tractation). On the other hand, the static threshold 
- computed with the same formula as that obtained in water - is well above 
the experimental data in the aeolian case: there is a very important hysteresis 
in that case. With both sets of data, we have thus a complete calibration of 
the dynamic and static threshold parameters, which allows to a good degree 
of confidence for the computation of the transport diagram on Mars (figure 4). 
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